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Abstract 

We consider slow motion of a pointlike topological defect (vortex) in the nonlinear Schrodinger equation minimally coupled to 
Chern-Simons gauge field and subject to external uniform magnetic field. It turns out that a formal expansion of fields in powers 
of defect velocity yields only the trivial static solution. To obtain a nontrivial solution one has to treat velocities and accelerations 
as being of the same order. We assume that acceleration is a linear form of velocity. The field equations linearized in velocity 
uniquely determine the linear relation. It turns out that the only nontrivial solution is the cyclotron motion of the vortex together 
with the whole condensate. This solution is a perturbative approximation to the center of mass motion known from the theory of 
magnetic translations. 
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Introduction 



I^^iis paper we are going to discuss the reliability of the Einstein-Infeld-Hoffman (EIH) method JjJ, which was developed in the 
tbferiry of general relativity for analysis of the motion of point sources of gravitational field, to dynamics of point-like topological 
delects in parity noninvariant systems. The essence of the EIH method, which remains relevant in the present context, can be 
outlined as follows. Let us assume that there is a characteristic velocity v in the model. In a relativistic model v is the velocity of 
light while in a nonrelativistic model v can be say a sound velocity. To consider time evolution which is slow as compared to the 
characteristic velocity, one can develope a perturbative expansion with X/v as an expansion parameter. The fields are expanded in 
powers of X/v and the time t is replaced by a rescaled time r, t = vt. The time derivatives are rescaled as d± = \d T . The n-th order 
time derivative of a given quantity is formally n orders of magnitude smaller than the quantity itself. In particular acceleration is 
always negligible as compared to velocity. However in a parity noninvariant system like a system in an external uniform magnetic 
field one can imagine a soliton performing a cyclotron kind of motion. In the cyclotron motion acceleration is always a linear form 
of velocity so that this kind of motion can not be described within the formal expansion in powers of X/v. 

The standard approach to the dynamics of point particles interacting with fields is to take in the first step arbitrary particle 
trajectories and find from field equations the fields produced by the assumed currents. In this first step the parameter X/v can 
be employed to calculate the fields an expansions in powers of particle velocities. The second step is to substitute the calculated 
fields into particles equations of motion. In this way one obtains purely mechanical equations which can be solved to obtain particle 
trajectories. Such a method was applied to derive the Darwin Lagrangian in electrodynamics Q or in the context of general relativity 
jlj to mention only the most important cases. 

It turns out that in the case of topological defects such an off-shell calculation is not possible. One can try to find solutions of 
field equations for a given soliton trajectory but it turns out that a regular solution exists but only for a very special trajectory 
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which is thus the solution of the soliton dynamics problem. This property has been observed first in the case of relativistic Chern- 
Simons vortices |J and independently for relativistic membranes [|J . The origin of the problem in the case of solitons can be traced 
back to the fact that, unlike for point particles interacting with fields, there is no a priori equation of motion for the topological 
defects. Unlike point particles solitons are composed of the same fields as the fields which mediate interactions between them. Their 
dynamics is implicit in the field equations. Our calculation in Section 3 provides another example. Its advantage is that there is an 
exact solution at hand to be compared with the perturbative result. 

2 Model and boost in an external magnetic field 

Let us consider the model ||, which is a field-theoretical description of the quantum Hall effect of polarized electrons. The electrons 
are described in terms of bosonic fields fermionized by coupling to the auxiliary Chern-Simons fields et M , 

L = ^ a %d a ap + J{**(A>*) - - -L(i? fc *)*p k *) - ^( Po - ***) 2 . (1) 

2 2 Ira 2 

Po, which is a condensate density, is related to the external magnetic field by po = —nB/e. In what follows we assume n to be 
negative. The covariant derivative couples the scalar field to both the external electromagnetic field and to the Chern-Simons 
field a M , = — iea^ — ieA^. The Greek indices run over space-time indices 0,1,2 while the Latin indices denote planar 
coordinates 1,2. Our convention is e 012 = +1 and we assume the signature (+, — , — ). In the formulation (|l|) the ground state of 
the theory in the external magnetic field B, 

A (t,x) = , 

A k (t,x) = -Be kl x l , (2) 

is the uniform condesate = p with the external magnetic field being screened by the Chern-Simons field, a k — — A k , &o = 0. 
Such a ground state admits topological vortex excitations || £| [Tj. 

The theory (0) is a Galilean invariant system in an external magnetic field. As such it has the following symmetry. If the set of 
fields ^>(t, x), a At (i, x) is a solution of the model (0), then the following boosted fields are also solutions of the model 

*(*,£) = *[t,x- R(t)]e lXB , 

a (t,x) = a [t,x-R(t)\ - R k a k [t,x - R(t)] , 

a k (t,x) = a k [t,x- R(t)] , 

1 . . . f* 

XB = --m{R k R k )t + mR k x k + e dr R k (r)A k [x - R{t)} . (3) 

2 Jto 

provided that the trajectory R(t) satisfies the equation of motion 

mR k = -eBs kl R l . (4) 

The last equation is simply the equation of motion of a planar electron in uniform magnetic field. Its solution is the cyclotron 
motion with the cyclotron frequency lo c = eB/m. If an unboosted solution is say the uniform condesate, then the solution after the 
boost is a condesate performing a cyclotron motion. If the unboosted solution contained a vortex then after the boost the vortex 
would move together with the condensate. 



3 Slow vortex motion in perturbative approximation 

In this section we are going to describe a perturbative scheme which can be applied in investigations of slow motions of topological 
defects in parity noninvariant systems. As the perturbative method is quite general, we think it may find applications in some 
other models. The method has been already described in our earlier article || devoted to the dynamics of relativistic self-dual 
Chern-Simons vortices. However, the calculation in || is relatively complicated and what is more there is no exact solution at hand 
to be compared with the perturbative result. The magnetic boost provides us with a nice example of exact solution, which confirms 
the perturbative result and the perturbative method as such. The perturbative calculation shows also uniqueness of the magnetic 
boost. 

Before we describe the perturbative calculation let us make a small rearrangement in the model (0). Namely, let us replace 
a-n + Ap = a M . This replacement and the use of the definition po = — — leads to the equivalent version of the model ([j]) 

L = ^ £ ^ ai d a ap - ep a + £{¥*(A,*) - *(^o*)*} - ^(r>**)*(U fc *) - \{ Po - ***) 2 , (5) 
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where we have already neglected the hats. The covariant derivative is simplified to = — iea^. The advantage of this 
formulation is that there is just the Chern-Simons gauge field to be handled with. The external magnetic field is replaced by the 
uniform background charge density now. 

We concentrate on the model @ from now on. The field equations of the model are 

id t ^ + ea * + ^-D k D k ^ + A(p - = , 

2m 

Ke k id k ai + e(p - p ) = , 

Ke kal3 d a a p + e J k = , (6) 
where p = is the particle density and the current is 

J k = ^-{*(D k *)* - = ^{d kX - ea k ) . (7) 

Ira m 

X is the phase of the scalar field, \P = yfpexpix ■ The model admits the uniform condensate solution ^> — ■Jpo, a M = — A^. Because 
there is a nonvanishing condensate, the model also admits topological vortex solutions. The Ansatz for a vortex solution with the 
winding number minus one can be taken as 

= \p{r)) l / 2 e- ie , 
af\S)=A(r) , 
a^^)=A (r) , 

4 0) (f) = . (8) 

The function p(r) interpolates between p(0) = and p(oo) = po- The gauge potential vanishes at the origin, A(0) — 0, and tends to 
a pure gauge at infinity, A(r) fts — In the Bogomol'nyi limit, A = jj^, the equations fulfilled by the profile functions in Eq.(||) 
can be derived || from 

2e 2 

V 2 lnp = (p- po) i 

K 

A(r) = -± + f , 
er Zep 

A (r) = -^-(p-po) . (9) 
Aran 

The primes denote derivatives with respect to r, ""' = Once the solution of the first equation in (||) is known the functions 
A(r) and Aq(t) can be expressed through p(r). 

The aim of the perturbative calculation is to find an approximate trajectory of a vortex in the limit of slow motion. The 
perturbative method consists of two main ingredients. The first of them is rather classic. As we are interested in slowly moving 
vortices, we expand all the quantities in powers of vortex velocity R k (t). The zero order approximation to a moving vortex solution 
would be just 

*(t,x) = * (0) [f-i?(i)] , 

a ll {t,x)=af[x-R{t)) . (fO) 

It is a solution to field equations but only when we neglect all the terms in field equations linear (or higher order) in vortex velocities. 
An exact solution $>(t, x), a Al (t, x), if exists, differs from the fields in Eqs.(|l0l), 

afj,(t,x) = a^(t,x) + u M (t,S) . (II) 

The deviations ip, are at least of first order in vortex velocity. The field equations (^), when linearized both in the vortex velocity 
and in the deviations ^,14^, become a set of inhomogenous differential equations 

idtt/j + ea ip + e*u + ^-[D k D k ^ - 2ie(D k V)u k - ie^(D k u k )} + Xp tp - 2A***V " A*V* = ~id t ^> , 
Ks k id k ui + e(* *V + * V) = , 



ka 



d a up -e 2 y*$Uk + e(d k x-ea k )(**ilj + -$ij*) = Ke M dta l , (12) 
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where = dk — ieak- The deviations if), u M can be expanded in powers of vortex velocity components R k (t). The leading terms in 
the expansion read 



ip{t,x) = R k (t)^[x- R(t)} 
u^{t,x)=R k (t)u^[x-R{t)] 



(13) 



Before we substitute the Ansatz (13) to the linearized equations (|12|), we have to introduce the second and key ingredient of 
the perturbative method. Namely we allow the vortex acceleration to be manifestly linear in velocity 



R k {t)=u kl R\t)+0{R k R k ) 



(14) 



The matrix to is time-independent and velocity- independent. In other words we are going to consider accelerations as being of the 
same order in magnitude as velocities. 

The assumption that the acceleration and the velocity are of the same order contradicts customary claims in the literature. The 
usual reasoning, which can be traced back to the Einstein-Infeld-Hoffman method in general relativity JO, is as follows. We assume 
there is a characteristic velocity in the model, say v. To consider slow time-evolution, we develope a perturbative expansion in the 
parameter l/v, which we consider to be small. The perturbative expansion has two ingredients. First of all the fields are expanded 
around boosted static solutions in the powers of l/v. Second, the real time t is replaced by a rescaled time r, i = vt. The time 
derivatives are then rescaled as dt = \d T . For a given quantity Q, its n-th time derivative is formally n orders smaller than the 
quantity itself, 9™Q = v~ n d™Q. In particular in this perturbative scheme the acceleration is always regarded to be negligible as 
compared to velocity or in other words the matrix u> kl in Eq.jT^) is implicitly assumed to be zero. Thus the commonly accepted 
perturbative scheme rules out any solution like a motion along circular orbit. 

We still do not have a satisfactory perturbative scheme to describe eventual vortex motion with radiation in which the back- 
reaction could be self-consistently taken into account. Thus in this paper we will restrict to nonradiative trajectories. When there 
is no radiation the energy of the vortex must be conserved. For a single vortex in an uniform condensate, translational invariance 
implies that its acceleration must be perpendicular to velocity 



to 



toe 



(15) 



where w is a constant. We do not know the value of the constant. It has to be fixed by a solvability condition. 
Let us then try to solve the Eqs.(|l^). The examination of sources on the RHS's of Eqs.(fl2|) 



id t V = — \R} cos 9 + R 2 sinO] 
2p 



-R 1 sin 9 + R 2 cos 01 



d t ax = [R 1 cos 9 + R 2 sin 01 A' sin 9 + [-R 1 sin 9 + R 2 cos 9} —cos 9 , 

r 

A 

d t a 2 = -[R 1 cos6 + R 2 sin 6] A' cos 6 - [-R 1 sin 6 + R 2 cos 9}- sin9 



(16) 



shows that without loss of generality we can adopt the following Ansatz for the first order field deviations (we assume the gauge in 
which the deviation of the scalar field's phase is zero) 



R k f^(x 
R k u { k) (x 
R k u {k) {x 
R k u { e k \x 



= p 1 1 ' 2 (r)s(r)[- R 1 sin 9 + R 2 cost 
= a(r) [-R 1 sin 9 + R 2 cos 9} , 
= bir^R 1 cos 9 + R 2 sin 9] , 
= c(r)[-7? 1 sin6» + 7? 2 cos6'] . 



(17) 



The substitution of the above Ansatz to Eqs.([l^) yields 



s" + 7 + —(s' - c) + 2a - 4ps = - , 

r r~ p r 

b' + + £■(&- - +2ws = , 

r r p r p 

c' + - - - - 2ps = , 

r r 

a — cub — p s — pc = — , 
r 

a 1 M 

luc — pb — A . 

r 



(18) 
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We use from now on the rescaled units in which all the constants m, (—k), pa and e are set equal to 1. We have also restricted to 
the Bogomol'nyi limit, which in the rescaled units corresponds to A = 1. It is by no means necessary but it makes the formulas 
more compact. In particular the second equation in the set (^), which holds only in the Bogomol'nyi limit, can serve for many 
simplifications. 

The second equation in the set ( |l8| ) is not independent, it can be derived from the last three equations. The last equation in 
the set ( fL8| ) can be used to express a(r) by other functions 

a = Lore + rpb + A'r . (19) 



Once a(r) is expressed like in Eq.(19), it can be eliminated from the first, third and fourth equation in the set ([18|). Finally we are 



left with only three independent equations 



b' + --- + 2us+^-[b-- + l]=0 , 
r r p r 

, c b 

c'+ 2ps = , 

r r 

s" + - - 4 + —(a' - c) - Aps + 2rpb + 2uorc = 2r{\ - A') . (20) 



p 



The regular asymptotes close to r = are 



s(r) w (1 + o)r + ... , 
b(r) ps a + ... , 

c(r) wa + ... , (21) 

where the higher order terms in r were neglected. There is only one free parameter a in the asymptote. 

To find out the asymptote at infinity we go back to Eqs.([l8|). At infinity p ~ 1, up to terms decaying exponentially. The 
asymptotic form of Eqs.(|T8|) is 

// s' s „ 2 

s" H j + 2a - 4s = - , 

, , b c 

b H h 2lus = , 

r r 

C ' + ^-2,s = , 
r r 

a' + — — tob + c , 

a 1 

This set of asymptotic equations implies the following 



s" + S - 
r 


s 


- 4s + 


n , a ' 
a H 

r 


a 


-2(1- 



luc . (22) 



(23) 

Once s(r) and a(r) are known, one can find b(r) and c(r) solving with respect to them the last two equations in the set (p2|). The 
solution is unique provided that lo 1 ^ 1. 

One can diagonalize Eqs.(^3|). The eigenvalues turn out to be 



A li2 = -2(l-w) . (24) 

If a; 2 > 1, one of the eigenvalues is positive and the other is negative. The asymptote must contain a linear combination of two Bessel 
functions and of two modified Bessel functions. Thus there is one divergent mode and two long-range unnormalisable modes to be 
removed. But we have only two adjustable constants, a in the asymptote at the origin (H) and lo. Thus there is no normalisable 
solution for u 2 > 1. 

\i lo 2 < 1, the general asymptote is a combination of two exponentially divergent modes and two normalisable modes. Thus the 
two constants u>, a might happen to be sufficient to remove the two divergent modes. Numerical analysis shows that it does not 
happen for — 1 < lo < 1, however. 
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Thus we must restrict our attention to the case uj 2 — 1. For to = +1 the asymptotic equations ( [22] ) are solved by 

a(r) w — + 2er , 
r 

e 2r g-2r 

s(r) « ^i— = +/3 2 — = +er , 
y2r v2r 

6(r) w -c(r) w 7 r~ 2 + |r 2 . (25) 

The two constants /3i and e must be tuned to zero with just one parameter a in the asymptote close to the origin. There is no 
regular solution for w = +1. 
For u> = — 1 the asymptote is 

, . 1 2a 

a(r) f» — I , 

r r 

e 2r „ e- 2r a 



s(r)«/3i -=+/?: 



V2r V2r r ' 
6(r) = -c(r) w 5 . (26) 

For nonzero <5 the gauge fields' deviations become R k af^ = SR l , compare with Eq.(p"7|). Once the constant /3i is tuned to zero with 
the help of the free parameter a, the asymptote of the electric current becomes J k = Se 2 poR k . Thus at any given moment of time 
the vortex velocity is parallel to the condensate current. What is more, for lo = — 1 one can give an explicit unique solution to 
Eqs.(pO|), namely b(r) = c(r) = —1 and s(r) = 0. According to Eq.(|l9|) 

a(r) =r[l - p(r) +A'(r)] = -A(r) . (27) 

The last equality holds thanks to the Gauss law. Thus the constant 5 is fixed as 6 = — 1. The vortex moves together with the 
negatively charged condensate. 

Note that u> = — 1 so the condesate moves along a cyclotron orbit, R k = —e kl R . The solution is a perturbative approximation 
to the exact deformed Galilean boost. Such a solution could not be obtained in a perturbative calculation if accelerations were 
neglected as compared to velocities. At the same time, the analysis shows that the deformed Galilean boost is the only solution 
which can be obtained by a perturbative expansion in velocities. In principle this does not exclude possibility of a dissipative 
noncyclotron motion. 
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